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We analyze the spin dynamics of an out-of-equilibrium large spin dipolar atomic Bose gas in an 
optical lattice. We observe a smooth crossover from a complex oscillatory behavior to an exponen¬ 
tial behavior throughout the Mott to superfluid transition. While both of these regimes are well 
described by our theoretical models, we provide data in the intermediate regime where dipolar inter¬ 
actions, contact interactions, and super-exchange mechanisms compete. In this strongly correlated 
regime, spin dynamics and transport are coupled, which challenges theoretical models for quantum 
magnetism. 

PACS numbers: 03.75.Mn , 67.10.Ba, 67.85.Fg, 05.70.Ln 


Dipolar atoms and molecules loaded in optical lattices 
are a promising platform to study quantum many-body 
physics HI, and in particular quantum magnetism [sj- 
[8|. In dipolar systems direct spin-spin interactions are 
provided by dipole-dipole interaction (DDI) without rely¬ 
ing on super-exchange mechanism |j3]. Although magne¬ 
tization changing collisions associated to the anisotropic 
character of dipolar interactions may introduce interest¬ 
ing exotic quantum phases 5M1 , these off-resonant 
processes are often negligible. Then dipolar interactions 
reduce to the following Hamiltonian: 
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where cP = /xo/ ^{gPBP (po being the magnetic per¬ 
meability of vacuum, g the Lande factor, gs the Bohr 
magneton), r is the distance between atoms, 0 i t 2 the an¬ 
gle between the magnetic field and the interatomic axis, 
and S±’ z are the spin operators acting on atom i. This 
Hamiltonian, known as the secular dipolar Hamiltonian 
in the context of Nuclear Magnetic Resonance 14j, bears 
strong similarities with the XXZ model of quantum mag¬ 
netism :S]. 

Experimental investigation of such spin Hamiltonians 
has recently started, with dipolar molecules 15|, mag¬ 
netic 


16] and Rydbergjl 
important interest 8, Lid 1 



atoms, which is triggering 
While these studies 


have focused in a localized regime where the particles are 
pinned to a well-defined position, in this paper we inves¬ 
tigate the case where magnetic atoms are free to move in 
an optical lattice. Thus spin dynamics and transport are 
coupled and conform to an intriguing interplay between 
super exchange mechanisms and dipolar spin exchange. 
Our experiment provides first data in this regime which 
challenges theoretical description. 

We study spin exchange dynamics of magnetic 


chromium 52 Cr bosonic atoms loaded in a 3D optical lat¬ 
tice, across the Mott to superfluid transition 0 . We 
observe, as a function of the lattice depth, a crossover 
between two distinct behaviors. In the Mott phase, spin 
dynamics displays a complex oscillatory behavior. We 
identify two distinct frequencies, one associated to on¬ 
site spin-dependent contact interactions and the other to 
inter-site DDI. In the superfluid regime, spin dynamics 
shows an exponential behavior which results from an in¬ 
terplay between contact and dipolar interactions. The 
regime at intermediate lattice depth is particularly in¬ 
teresting because super-exchange mechanisms also con¬ 
tribute to the spin dynamics. For example an atom may 
tunnel into an already occupied site and interact with 
another atom by spin dependent contact interaction; this 
can trigger spin changing collisions (see Fig(T}^)) . It is ex¬ 
tremely challenging to simulate the many-body quantum 
spin dynamics in this intermediate regime where three 
exchange mechanisms compete (associated to DDI, con¬ 
tact interactions, and super-exchange). We experimen¬ 
tally find that the oscillations observed in the Mott phase 
survive with a reduced amplitude. 


We perform our experiment with a spin-3 52 Cr Bose- 
Einstein condensate, compri sing 10 4 atoms loaded into 
an anisotropic optical lattice 16]. As the lattice depth is 
spanned, we observe the superfluid to Mott transition, at 
a typical lattice depth of 12 Er (where Er is the recoil 
energy). For our experimental parameters, the system 
in the Mott phase consists of a core with two atoms per 
site, surrounded by a shell with one atom per site. 


To initiate spin-dynamics atoms are transfered into the 
first single-particle excited Zeeman state m s = — 2 using 
the tensor light-shift of a 427.85 nm light pulse [I(|. We 
then measure, after a variable hold time t, the spin popu¬ 
lations by means of a Stern-Gerlach procedure. A typical 
result is plotted in FigHJ:). Only the m s = —3, —2, — 1 
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FIG. 1: a) Simple representation of the system in the Mott phase with one atom per site, with atoms interacting via DDI 
(white ellipse); b) sketch illustrating competition between exchange due to DDI (Vdd) and tunneling (J) assisted spin-exchange 
due to contact interactions ( 14 ); c) typical measurement of the spin components as a function of time; d) time evolution of 
observable n_ 3 /n _2 for two extreme regimes (top : superfluid, bottom : Mott). Lines are guides for the eye resulting from fits. 


and 0 spin components are significantly populated as the 
system evolves. The populations display a rather com¬ 
plex behavior as a function of time. In order to simplify 
the discussion, we focus our attention onto the observable 
given by the ratio n_ 3 /n _2 of m s = —3 and rn s = —2 
populations, since they are the most populated compo¬ 
nents. 

We have plotted in Fig[l] d) the typical results cor¬ 
responding to the extreme regimes of high and shallow 
lattice depths, showing quite different spin dynamics. In 
the Mott phase we observe at short times (< 0.5ms, see 
inset) a strongly damped oscillation, and then at longer 
times a second oscillation. In the superfluid phase the 
spin dynamics is better described by an exponential. All 
these features are present in the data from 3 Er to 25 
Er. We plot in FigEJs, the amplitudes of the exponential 
and of the slow oscillation, and in FigEJo the frequencies 
of the fast and slow oscillations. When reducing the lat¬ 
tice depth, we observe that both oscillation frequencies 
decrease and they become closer to each other. Oscilla¬ 
tions do survive at lattice depth slightly below the Mott 
to superfluid transition. However, for very shallow lattice 
depth, the oscillations at low frequency are barely visible 
and the spin dynamics is mostly exponential. 

For most of the data shown in Fig. El magnetization 
is constant (see the inset of FigEO). The stability of 
the magnetization for large lattice depths indicates that 
dipolar relaxation is suppressed, which arises because the 
energy released in a dipolar relaxation event (the Larmor 
energy rs 30 kHz) does not match band excitation (> 
50kHz) [2(J. At small lattice depth, dipolar relaxation 
is not completely suppressed (as the first excited band 


has an energy close to the Larmor energy). However, 
this effect remains rather small, and is neglected in the 
theoretical analysis presented below. 

To account for the dynamics at the lowest lattice 
depths, the Gross-Pitaevskii equation can be safely used 
because the gas is in the condensate phase. We performed 
a numerical simulation up to a lattice depth of 7 Er in or¬ 
der to describe the observed data. The interaction term 
takes into account short range contact interaction and 
non local DDI 21]. Concerning DDI we only include the 
spin conserving terms (see Eq. ©)• 

Simulations display a complex behavior, see Fig. [3] for 
7 Er. The general trend of the spin dynamics, showing a 
slow drift with a s=s 14 ms characteristic 1/e time, is well 
reproduced by the simulation. By fitting both curves 
(experimental and numerical) by an exponential, a good 
agreement for the amplitude is found (see Fig.Ek))- 

In addition, Fourier analysis of the numerical results 
(up to 15 ms) displays a lowest resolved frequency in good 
agreement with the largest experimentally observed fre¬ 
quency at a lattice depth of 3 Er and 7 Er (see Fig.Eb))- 
Nevertheless, the strong damping observed in the exper¬ 
iment is not reproduced by the zero-temperature simu¬ 
lation. For very low lattice depths (below 3 Er), simu¬ 
lations show that both contact and dipole interactions 
contribute to the spin-dynamics. In particular if we set 
DDI to zero we numerically observe a spin dynamics fre¬ 
quency which is roughly twice faster. This illustrates the 
interplay between dipolar and contact interactions in the 
superfluid regime. 

We also investigate numerically the spatial dependence 
of the dynamics in Fig. Eh>) , which shows a cut of the den- 
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Lattice depth (Er) 

FIG. 2: Spin dynamics amplitudes and frequencies as a func¬ 
tion of the lattice depth, a) Amplitude of the exponential dy¬ 
namics (blue diamonds) and slow oscillation (red triangles). 
Green circles are results of numerical simulations. Inset: vari¬ 
ation of the magnetization over 20 ms. Solid lines are guides 
to the eye. b) Frequency of fast (black points) and slow (red 
triangles) oscillations. The black top solid line corresponds to 
spin-exchange frequency associated to intrasite contact inter¬ 
actions, while the black open circles correspond to a numeri¬ 
cal simulation of Gross-Pitaevskii equation. The bottom red 
curve is a guide to the eye. The blue dot-dashed line shows 
the prediction in the Mott regime (see text). The frequency 
of super-exchange process is given by the green solid line. Er¬ 
ror bars in frequency and amplitude result from the statistical 
uncertainty in the fits. 


sity of the atoms in the m s = —2 Zeeman state in absence 
of an optical lattice, for three evolution times. Similar re¬ 
sults are obtained at small lattice depths. The dynamics 
is inhomogeneous, due to a nontrivial interplay between 
contact and dipolar interactions. While spin-exchange 
interactions due to contact interactions are larger in the 
center of the cloud, due to higher density, dipolar in¬ 
teractions are stronger in the outskirts 22]. The ob¬ 
served dynamics is faster in the center, which illustrates 
the dominant role of contact interactions at low lattice 
depths. 


We now turn to our theoretical analysis at large lat¬ 
tice depth, where the system is not superfluid and one 
should go beyond the Gross-Pitaevskii mean-field de¬ 
scription. We focus on the frequencies of the two os¬ 
cillations observed in the Mott phase. We stress that 


n,/n 2 (a) 



Time (ms) 



FIG. 3: Results of numerical simulation using the Gross- 
Pitaevskii equation, a) Spin dynamics for the deepest lattice 
depth that could be simulated (7 Er). The (red) triangles 
are experimental data, b) Spatial analysis of spin dynamics, 
showing a cut of the density of the atoms in the m s = —2 
state (with no lattice on). 


for the largest lattice depth, the experiment is performed 
in a regime where tunneling is practically absent, and 
super-exchange interactions are exponentially reduced, 
as shown in Figj2}o) . 

We interpret the upper frequency as a result of intra¬ 
site spin-exchange dynamics (| — 2 , — 2 ) —I ^ (| — 3, — 1 ) + 
| — 1,-3)) arising from spin-dependent contact interac¬ 
tion in doubly-occupied sites. The observed frequency 
is in good agreement with the theoretical frequency 


2 V 2 m — ° 4 ) n o ( see black solid line in Figl2b))- Here, 
no is the peak density in a doubly-occupied lattice site, m 
the mass of the atoms; ciq and a 4 the scattering length of 
S = 6 and 5 = 4 molecular channels respectively 23[ 24]. 
As also shown in Figl^bO, weaker on-site confinement in 
shallower lattices reduces the density n 0 ; and therefore 
the observed frequency. 

As for the lower frequency, we associate it to non-local 
DDI between doublons. Indeed, as reported in our pre¬ 
vious work jl 6 j |. this oscillation disappears in presence 
of a strong magnetic field gradient, which proves that 
it originates from inter-site interactions; in addition, the 
oscillation is also absent when doubly occupied sites are 
emptied [ 211 . 

We will now describe a perturbative model to account 
for the many-body interactions between doublons due to 
intersite DDI. We first assume that, after the damping of 
the fast oscillations, the ensemble of doublons is in a sta¬ 
tistical mixture of |S = 6 ,m s = 4) and |S = 4, m s = 4). 
We then describe the spin-dynamics due to inter-site DDI 
between doubly occupied sites by the following model, 
inspired from Gm. We calculate the time evolution 
of the population 7V _2 in the state m s = —2 using per¬ 
turbation theory in the Heisenberg picture. The many- 
body Hamiltonian takes into account the interaction of 
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one doubly occupied lattice-site i with all its neighbors j 
by pair-wise DDL Using Taylor expansion, the formal ex¬ 
pression of the population reads N_ 2 {t) = M n t n . 

While Mi and M 3 vanish, the second moment M 2 de¬ 
scribing spin dynamics up to the second order of pertur¬ 
bation reads: 


M - 2 = -Y J v dd{n, J ) 2 /h 2 


( 2 ) 




where Vdd(u,j) is the dipolar spin-exchange matrix ele¬ 
ment between sites i and j [ 21 ]. 

From M 2 , we first extract an estimate of the spin os¬ 
cillation frequency v. cos 2 (nist) « 1 — M 2 f 2 . However, 
as shown in 14], the higher order terms are expected to 
lead to a reduction of spin dynamics amplitude and a de¬ 
crease of the quasi period by a factor of typically two. We 
indeed recover these features by taking into account the 
moment up to the fourth order M 4 (following the calcu¬ 
lation in 14j). We apply this perturbative approach up 
to fourth order to the case of an assembly of doublons in 
state 15 = 6 ,m s = 4), for which intersite spin-exchange 
interactions are strongest. 

The result is displayed in Fig On), and shows good 
agreement with the observed frequencies, which illus¬ 
trates the relevance of our model in the deep Mott regime. 
The observed reduction of the frequency at lower lat¬ 
tice depth may be the consequence of a reduction of 
the doubly-occupied sites Mott plateau, leading to either 
stronger border effects, or more defects. The reduction of 
the frequency in the presence of holes is also an expected 
feature 0 - The damping (also observed in 0 ) is not 
reproduced by our theoretical model, which can only ac¬ 
count for dynamics at short times. In general, it should 
be stressed that the many-body system which we study 
here is extremely challenging for non-perturbative many- 
body simulations, due to the large spin and the immense 
Hilbert space which has to be taken into account. 

In this paper we have explored spin dynamics of 
chromium atoms loaded in an optical lattice as a function 
of the lattice depth. Superfluid and Mott regimes lead to 
markedly different features. Experimentally, the spin dy¬ 
namics evolves smoothly when we vary the lattice depth. 
However, our analysis shows that the impact of DDI on 
spin dynamics is drastically different in the two regimes: 
whereas in the weakly-interacting regime, dipolar inter¬ 
actions are described by a mean-field associated to a 
geometrical average, in the strongly correlated regime, 
dipole-dipole couplings from one site to the other sites 
add quadratically. One important consequence is that 
spin dynamics mediated by dipolar interactions should 
be a border effect in the mean-field regime, but a bulk 
effect in the strongly correlated gas 21]. Space-resolved 
measurements could therefore be a very interesting way 
to discriminate both regimes. 


Whereas the two extreme regimes of shallow and deep 
lattice depths are qualitatively well described by our the¬ 
oretical models, one of the novelties of the paper is the 
experimental study of the regime close to the Mott to 
superfluid transition. In this case, exchange processes 
due to dipolar interactions, spin-exchange due to con¬ 
tact interactions, and super-exchange interactions, may 
all contribute to the dynamics on the same footing. Our 
experiment provides first insights into the coupled out-of- 
equilibium magnetic and transport properties of such a 
strongly correlated gas, which challenges theoretical de¬ 
scriptions. Our study thus pioneers the study of the mag¬ 
netic properties of an exotic novel quantum many-body 
system made of large spin dipolar particles. 
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Supplemental material 


GPE Simulation 

To account for the dynamics at the lowest lattice 
depths, we use the Gross-Pitaevskii equation. We per¬ 
formed a numerical simulation in order to interpret the 
observed data. The Gross-Pitaevskii equation for this 
system reads: 

ih^ n (t,r) = (H° tl + 

+ Jp q ,k{t,r')V ntqtkt i(r-r')d 3 r , ')'ipi(t,r) (3) 

where ip n (t, r) is the wave function for the nth spin com¬ 
ponent in the z direction, p q j~{t,r) = is 

the one-body density matrix and V n ^k,i(r) is the interac¬ 
tion term. The summation on repeated indices is implied. 

The term if ° ; is the single particle term of the Hamil¬ 
tonian and reads 

Ki = S n ,i + P ext (rl) + gpsBS^ (4) 

This term includes the kinetic energy and potential en¬ 
ergy which contains the external trapping and the mag¬ 
netic potentials. In order to take into account the lattice 
potential we renormalize the mass [J . 

The interaction term takes into account short range 
contact interaction and DDL The short range one 
is written as I4 r = Ej =0 2 4 69i^i> where Pj = 

Yln=-j b'i n )0i n \the projector on the molecular state 
of total spin j. To account for the lattice, we use 
9j = ( 27 r)V 2 fj 3 . where gj = aj (aj being the scatter¬ 
ing length in channel J), and a is the geometrical average 
width of the lattice site. 

Concerning DDI we only include the spin conserving 
terms since in presence of a lattice the dipolar relaxation 
is strongly inhibited [2] (see Eq. ©)• Dipolar interac¬ 
tions have a non-local character which is explicitly taken 
into account in Eq. (J3|). Due to geometrical averaging, 
dipolar interactions are typically smaller in the center of 
the cloud compared to the outskirts. Dipolar coupling 
constant does not need to be renormalized. 

Super-exchange 

The super-exchange mechanisms in the Mott regime 
for a spinor gas become rather involved due to the multi¬ 
channel character of the interaction. A large number 
of additional spin-exchange phenomena can thus arise in 
addition to the exchange of the spins found in spin 1/2 
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Estimated doublon fraction 


Spin-exchange as a function of temperature 


To demonstrate that the slow spin oscillations which 
are observed in the experiment at large lattice depth are 
due to interactions between doubly-occupied sites, we 
have varied the number of doublons by modifying the 
temperature of the gas before it is loaded in the optical 
lattice. This is achieved by tuning the value of the trap 
depth at the end of forced evaporation. While the tem¬ 
perature is always kept below the energy gap for band ex¬ 
citation in the optical lattice (corresponding to typically 
2.5 ftK), the initial temperature has a strong impact on 
the atom number distribution in the lattice sites. 



0 5 10 15 

Time (ms) 


FIG. 4: (a) Oscillation amplitude as a function of the fraction 
of doublons in the cloud. The data point at the origin corre¬ 
sponds to a situation where all doublons are removed using 
dipolar relaxation at large magnetic field. The other points 
correspond to different temperatures, (b) For each point, the 
doublon fraction is independently estimated by measuring the 
fraction of atoms which disappear when doubly occupied sites 
are emptied using dipolar relaxation at large magnetic fields 


systems. For example, for a spin 3, a two-body state |2,1) 
(where one site contains an atom in state with m s = 
2, and the next site contains an atom with m s = 1) 
could become 11, 2) or |3, 0) or |0,3). The matrix elements 
describing tunneling-assisted super-exchange with spin- 
dependent contact interactions read: 


T T m' i ,m ' 0 
U m \,7712 


E 

S,ms 


-—(m 1 ,m 2 \S,m s )(S,m s \m , 1 ,m' 2 ) (5) 
Us 


where \S,ms) is the molecular state of total spin S with 
projection ms, and (S, ms\m[, m' 2 ) are Clebsch-Gordan 
coefficients. mwi t2 ) represent atomic magnetic states. 
J is the (mi independent) tunneling matrix element be¬ 
tween nearest sites, and Us the spin-dependent on-site 
contact energy (which only depends on S). 


For each given temperature, we characterize the lattice 
atom distribution by measuring the fraction of atoms lo¬ 
cated in doubly-occupied sites. To this end, we use the 
following protocol. The magnetic field is set to a value 
> 1G such that the Zeeman energy significantly exceeds 
the trap depth. Atoms are promoted to the Zeeman ex¬ 
cited state m s = 3 using a 5 ms radio-frequency sweep. 
At such a large magnetic field, dipolar relaxation between 
a pair of atoms is a local process Q, which releases an 
energy that is large enough for both atoms to leave the 
trap. Then, as shown in Fig. S](b), the number of atoms 
in the lattice rapidly decreases, towards a steady-state 
non-zero value. This is due to dipolar relaxation occur- 
ing in doubly occupied sites, while isolated atoms are 
immune from dipolar relaxation. The measurement of 
the fractional loss after dipolar relaxation is therefore a 
measurement of the fraction of atoms located in doubly 
occupied sites. 

For each temperature, we have also performed a spin- 
dynamics experiment identical to the one described in the 
main part of the paper (i.e. after atoms are promoted to 
state ms = —2). We then observe a similar slow spin 
oscillation, whose amplitude is plotted in Fig. |4ji) as a 
function of the estimated number of doublons. This fig¬ 
ure demonstrates that the oscillation which is seen in the 
experiment has an amplitude which increases with the 
number of doublons. As also reported in Q, we observe 
no oscillation when the doubly occupied sites are emp¬ 
tied using dipolar relaxation at large magnetic field be¬ 
fore studying spin dynamics from state m s = —2. Taken 
together, these measurements therefore confirm that the 
slow oscillations are due to the interactions between dou¬ 
bly occupied sites (the non-local character of the inter¬ 
action at play is demonstrated in !)• 

We also stress that the measurements shown in Fig [4] 
(a) indicate that spin oscillations survive at relatively 
large temperature, corresponding to clouds above the 
BEG critical temperature before they are loaded in the 
lattice. The oscillations which we see in our experiment 
thus survive despite the presence of the disorder in the 
number of atoms in different lattice sites which necessar¬ 
ily occurs at non-zero temperature. 
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n,/n 2 



Time (ms) 


FIG. 5: Time evolution of observable n_ 3 /n _2 for a 2D array 
of ID quantum gases produced by loading atoms in a deep 
2D optical lattice. 


Spin-exchange in a 2D lattice 

The main part of this paper presents the study of spin- 
exchange dynamics due to intersite dipolar interactions 
for a range of lattice depths spanning the Mott to super- 
fluid transition. We find that the situation is particularly 
original and interesting in the superfluid regime, as spin 
dynamics is then driven by an Heisenberg-like Hamilto¬ 
nian despite the fact that the atoms are not in an insu¬ 
lating state. Unfortunately, we find that magnetization¬ 
changing collisions due to dipolar relaxation are not com¬ 
pletely negligible at the lowest lattice depths (see the in¬ 
set of Fig©.)). This is because the Zeeman energy is then 
similar to the gap for band excitation, so that dipolar re¬ 
laxation is energetically allowed [2j. Although the agree¬ 
ment between our experimental data and the mean-field 
simulation based on the Gross-Pitaevskii equation ne¬ 
glecting dipolar relaxation is still satisfactory, the parallel 
to the Heisenberg model of magnetism may then be ques¬ 
tionable, because the Heisenberg model of magnetism 
does not include magnetization-changing collisions. 

To probe spin-dynamics in a mean-field (superfluid) 
regime, while maintaining a constant magnetization in 
time, we have therefore performed additional similar ex¬ 
periments in a deep 2D optical lattice, where dipolar re¬ 
laxation is strongly reduced due to a small density of 
states at the Larrnor energy [s], Then the atoms form 
a periodic array of one-dimensional quantum gases. As 
shown in figure 0 after the atoms are promoted to state 
ms = —2, we observe spin dynamics with an oscillatory 
behavior, at constant magnetization, similar to what is 
obtained in a 3D optical lattice. The oscillation frequency 
which is deduced from a fit using a damped sinusoidal 
function, 165 Hz, is significantly smaller than the fre¬ 
quency measured in the deep 3D optical lattice, but com¬ 
parable to the ones measured for lower 3D lattice depths. 

The theoretical interpretation of spin-dynamics in this 
2D lattice is difficult, because dipolar interaction may in¬ 
clude intersite couplings between the ID quantum gases 


created by the optical lattice, and interactions within one 
tube of atoms. In addition, spin-dynamics due to on-site 
spin-dependent contact interactions may not be negligi¬ 
ble. However, the fact that the observations are quali¬ 
tatively the same above the 3D Mott transition in a 3D 
lattice and in an array of ID gases in a 2D lattice is an 
interesting feature, which shows that the basic out-of¬ 
equilibrium magnetic properties of the insulating state 
are preserved in presence of transport. 


Dipolar interactions in the mean-field and in the 
strongly correlated regimes 


As shown in the main part of the paper, dipolar in¬ 
teractions lead in principle to drastically different be¬ 
haviors in the mean-field and in the strongly correlated 
regimes. In the mean-field regime, dipolar interactions 
between atoms create a non-local mean-field, determined 
by a global average over space 

4>(r) oc J d 3 r'V d d{r - r')n(r') (6) 

To simplify this specific discussion, the tensorial nature 
of Vdd is overlooked; see Eq. ©. Eq. © in practice 
implies that for a homogeneous system, dipolar interac¬ 
tions vanish and should not contribute to spin dynamics. 
Therefore, spin-dynamics due to dipolar interaction is a 
border effect in the mean-field regime. 

On the other hand, in the Mott regime, according to 
second order perturbation theory applied to an initial 
state where each lattice site is populated by one atom in 
state m s = —2, spin dynamics due to dipolar interactions 
is qualitatively described by the following effective rate 
at short times 


r = ^/E^b) 2 (7) 

V *<j 

where 

d 2 (l -3 cos 2 (6^) 

VddKTxj) — 4 3 

hj 

<3, —2; 3, —2| S+ S~ |3, —1; 3, —3) (8) 

describes spin exchange dipole-dipole interactions be¬ 
tween two m s = —2 atoms in sites i and j. When we 
analyze the case of doubly occupied sites, we consider 
spin-exchange between two S = 6 pairs of atoms in state 
IS = 6, ms = —4), and we use 


VMn, t ) = x 


4r? • 

i,3 


(6, —4; 6, —4| Sf S~ 16, —3; 6, —5) 


( 9 ) 


The difference between Eq. [6] and Eq. [7] is striking and 
illustrates the difference between the two regimes. In the 









superfluid regime, the effect of dipolar interactions results 
from a geometrical averaging, whereas in the strongly 
correlated regime, it results from a quadratic averaging. 
As a consequence, in the strongly correlated regime, spin- 
exchange due to dipolar interactions does not vanish even 
for a homogeneous system. 

To better understand the difference between Eq. [G]and 
Eq. 0 we propose the following physical interpretation. 
The initial many-body state is ’I'o = |— 2, — 2, — 2), 
where the spin state of each atom in each lattice site is 
explicitly written down: each site contains one m s = —2 
atom. Dipolar interaction is a sum of pair-wise interac¬ 
tions (see Eq. 0 between atoms in sites i and j. There¬ 
fore is directly coupled to the following many-body 
state: 

'JT = , 1 = x 

Vd,d{ r ij ) 2 

Y.Vddinj) |-2,..., i : -1 ,..,j : -3...,-2) (10) 

i<j 

In the state |—2: —1 ,..,j : —3...,—2) all sites con¬ 
tain a m s = —2 atom, except site i (resp. j ) which 
contains a m s = — 1 (resp. -3) atom. The rate of cou¬ 
pling between 'To and is given by Eq. 0 This many- 
body physical picture clearly demonstrates that the ini¬ 


tial state is coupled to a state which shows spin entan¬ 
glement. As a consequence, spin dynamics in this regime 
cannot be grasped by the Gross-Pitaevskii equation and 
by Eq. EO 

This discussion indicates that the difference between 
the regime at low lattice depth and the situation at 
large lattice depth could be further experimentally in¬ 
vestigated, either by measuring the growth of entangle¬ 
ment (for large lattice depth), or by performing position 
resolved spin dynamics measurements. 
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